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Investigations of the coordinates in Ptolemy’s 
Geographike Hyphegesis Book 8 


Christian Marx* 


Abstract In Book 8 of his Geographike Hyphegesis Ptolemy gives coordinates for ca. 
360 so-called noteworthy cities. These coordinates are the time difference to Alexan¬ 
dria, the length of the longest day, and partly the ecliptic distance from the summer 
solstice. The supposable original conversions between the coordinates in Book 8 and 
the geographical coordinates in the location catalogue of Books 2-7 including the un¬ 
derlying parameters and tabulations are here reconstructed. The results document the 
differences between the 17- and S-recension. The known difference in the longitude of 
Alexandria underlying the conversion of the longitudes is examined more closely. For 
the ecliptic distances from the summer solstice of the fl-recension it is revealed that 
they were originally computed by means of a so far undiscovered approximate, linear 
conversion. Further it is shown that the lengths of the longest day could be based on 
a linear interpolation of the data in the Mathematike Syntaxis 2.6. 

Keywords Ptolemaios, Geographike Hyphegesis Book 8, Conversion of coordinates, 
Longitude, Latitude, Length of the longest day. Ecliptic distance 


1 Introduction 


Ptolemy’s Geographike Hyphegesis (GH/Geography) compiled in about 150 AD con¬ 
tains a location catalogue in Books 2-7 with about 6,300 ancient localities and their 
positions expressed by geographical longitude A and latitude <&. 360 of these localities 
are additionally listed with coordinates in Book 8 of the Geography. In contrast to the 
location catalogue, the positions in Book 8 are expressed by means of 

- the time difference A (in hours) from the location to Alexandria, 

- the length of the longest day M (in hours) at the location, 

- the ecliptic distance S (in arc degree) of the sun from the summer solstice at the 
time when the sun reaches the zenith (for locations between the tropics). 

The localities in Book 8 are selected cities, the so-called poleis episemoi (“noteworthy 
cities”). They are arranged in 26 chapters; each chapter is related to one of the maps 
Ptolemy divided the Oikoumene into (the inhabited world known to the Greeks and 
Romans). Figure [J shows a gene ral view of the sites of the single maps (see also 


Stiickelberger and Grafih^ 2006ll . 


Since the positions of the poleis episemoi are given in the location catalogue and in 
Book 8, a comparison of the positions of both sources is possible after a conversion 
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of coordinates. The poleis episemoi also appear in another of Ptolemy’s works, the 
Procheiroi Kanones (Handy Tables; hereinafter denoted Canon). Unlike Book 8 the 
coordinates are given there in form of geographical coordinates as in the location 
catalogue of the Geography. 

A convers ion of t he A- and M -data in Book 8 into geographical coo rdinates can 
be found in Ciintd ( 192,‘lh and in [Stiickelberger and Mittenhubeij ( 2009ll . The latter 
contains an edition of the Canon with a compilation of the corresponding geographical 
coordinates of the GH 2-7 and the converted values A(A) and $(M) on the basis 
of Book 8. Cunt j examines four maps of Europe (GH 8.5-8). He supposes that A 
and M in Book 8 originate from the geographical coordinates in Books 2-7. By his 
assumptions about the precision of A and M as well as about tabulations underlying 
A and M he makes an atte mpt to explain the coordinate values in Book 8 (see Sects. H] 
and ED- iHonigmannI (|l929l . p. 64ff.) assumes that the data of the poleis episemoi in 
Book 8 as well as in Books 2-7 were derived from works by a predecessor of Ptolemy. 

In the following, the conversions of coordinates between the location catalogue and 
Book 8 are investigated. The underlying calculation methods and their parameters are 
worked out as well as the underlying tabulations. Deviations between the coordinates 
in Book 8 and in the location catalogue are explained as comprehensively as possible by 
the resolutio n (precision) of the coordinates. As source for the data of the Geography, 
the edition of Stiickelberger and Grafiho^ ( 2006ll was used (for changes on the data see 
Appendix). It contains the coordinates of the two main recensions of the Geography, the 
U- and S-recension. The latter one has been handed down only by the X-manuscript; 
its location catalogue ends in Book 5, Chapter 13, so that a direct comparison with 
the data in Book 8 is not possible . Both U and X are considered in the present work. 
Stiickelberger and Grafihofi ( 20061 . p. 44) state that many M-values of X seem to be 
based on the coordinates of the Canon. The relation between the coordinates in Book 
8 and the Canon is only broached here. The coordinate s used are based on the edition 
of the Canon bv lstiickelberger and Mittenhubei (200^. 


2 Time difference to Alexandria 


The conversion between the time difference to Alexandria A (in hours) and the longi¬ 
tude A (in degree) is 

A = A 15Vh-HAA (1) 

with the longitude of Alexandria Aa (in degree). Ptolemy gives two different values 
for Aa in the Geography, in the location catalogue 60°30' =: Aa (4.5.9) and in Book 
8 4h =: Aa (8.15.10; the distance between Alexandria and the zero meridian); in 
GH 7.5.14, he gives both values. The difference between both parameters is AAa := 
A a - Aa = 60° 30' - 60° = 30' = 2 min. 

ICunt j ( 19231 p. 101 ff.) investigates the coordinates of the 3rd to 6th map in GH 8.5- 
8 (Europe, 42 localities). He states that both 60°30' and 60°00' underlie the time 
differences Ai {i = 1,2,... ,number of coordinates n) and that it is not possible to 
determine the Aa the coordinates are based on in all cases. His explanation for different 
Aa is that Ptolemy upgraded the value 60° 30' to 60° 00' and added the updated Ai in 
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a furth er column besides the values under lying 60° 30'; later copyists chose one of the 
values, fstiickelberger and Grafihofil ( 20061 p. 44) mention that 60° primarily underlies 
the data of as well as 60°30' the data of X. 


2.1 Conversion of coordinates and accuracy 

Ptolemy explicitly declares in GH 8.2.1 that the time differences Ai in Book 8 are 
approximate values; thus, equality after a conversion between the data of Books 2-7 
and Book 8 is not to be expected. For the following conversions of coordinates, the 
difference Ay between the conversion y{x) of a coordinate x and the given coordinate 
value y is defined as 

Ay = y - y{x) . (2) 

Then Ay is the error of y with respect to y{x). 

Both Aa and Aa were used here for the conversion of coordinates. For the conversion 
A !->• A by Eq. ([T]) Fig. [5] shows the frequencies of differences AA^ = — A(Ai) lying 

within specified intervals (interval width 10'). For fl the conversion with Aa gives a 
maximal frequency near zero and with Aa near —30'. That is the negative systematic 
difference AAa and indicates that the majority of the coordinates were originally 
converted with Aa = 60°. In X the conditions are converse so that Aa = 60° 30' 
must be the main parameter of the original conversion. The frequency distribution is 
bimodal; the distance of the lower second peak from the main peak is 30'. Hence, the 
second peak is probably caused by Ai which are based on the corresponding other Aa. 

The size of the differences AA^ can mainly be explained by a original conversion 
A I—>■ A in conjunction with a rounding of A. Apart from a few excepti o ns, th e largest 
denominator of the fractions of the Ai is 30 (1/30 h = 2 min). ICunt j ( 19231 p. 108) 
assumes an original precision of 1/12 h = 5min (= 1°15')0 In general, if a number is 
rounded with a precision (resolution) a = l/d, the expected maximal rounding error 
is ^ 

emax = ±- = ±— . (3) 

For a = 5 min Cmax is 2.5 min 2 / 3 °. This agrees with that part of the tail of the 
frequency distribution which is not extended by the systematic difference AAa (see 

Ic^t j I I 923 L p. 108) gives a tabulation of A- and A-values with a step size of 


AA = 5min (AA = lV4°) by which A-values may have been determined. ICuntz 
supposes that in the case of a A^ which does not equal a tabulated A the nearest 
tabulated A was chosen. In the end I Cunt a cann ot expla in several Ai by means of the 
tabulation. The use of a tabulation described bv I Cunt 3 was applied here to the entire 
A-data of Book 8. A tabulation was generated for Aa and Aa. In the case of H 61% 
and in the case of X 51% of the A-values determined by one of both tabulations equal 
the Ai given in Book 8. Furthermore a tabulation with A A = 2.5 min wa s tested ; the 
result are 61% and 55%. That does not argue for the use of a tabulation. Ciintd gives 
the example of Lugdunum (8.5.5) with A = 23°15' and A = —2’°30™. The nearest 


iHonigmannI lll929l '). who discusses ICunt j ’ work, does not dissent. 
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tabulated A is 23° so that A becomes —2*'30™. However, this might also be the result 
of Eq. (P) and a rounding: A(A) = -2^29“. 

2.2 Determination of Aa and the resolution 


Considering the single chapters of Book 8, there are chapters whose longitudes are 
converted with one of both Aa mainly; in other cases, the underlying Aa cannot be 
recognized reliably. The question of the Aa underlying a single A led to an attempt to 
determine Aa for each locality by means of the following method. A^ is converted with 
Aa to A(Ai) and with Aa to A(Ai) by Eq. (p. Ai and Ai are rounded with different 
precisions a to determine the apparent resolution of the given Ai. In general terms, the 
two coordinates x and y of different types are given as well as the converted coordinate 
y{x). Eurthermore a main resolution Og is set. The procedure applied is then 
Procedure ROUND: 

(1) Round y{x) to the precision as] if Ay = y — y{x) = 0 end. 

(2) For each of the precisions a < as downwards: round y(x); if Ay = 0 end. 

(3) For each of the precisions a > Og upwards: round y(x); if Ay = 0 end. 

If Ay becomes 0, the apparent resolution of y is found; if Ay does not become 0 until 
the end, an erroneous y ot x can be assumed. 

Procedure ROUND was applied to the Ai and Ai. The results of a comparison with 
the given Ai were valued as follows: 

(1) rounded Ai = Ai or/and rounded Ai = Ai: 

(a) only rounded Ai = Ay. Ai is apparently based on Aa; 

(b) only rounded Ai = Aii Ai is apparently based on Aa; 

(c) rounded Ai = Ai and rounded Ai = Ay. AA^ = — A^, AA^ = A^ — A^ 

(i) IA Ail < |AAi|: Ai is possibly based on Aa; 

(ii) IA Ail > |AAi|: Ai is possibly based on Aa; 

(in) IA Ail = |AAi|: no decision possible; 

(2) rounded Ai A Ai and rounded Ai A Ai: no decision possible, maybe a corrupted 
Ai-value. 

The possible resolutions for Procedure ROUND were chosen according to the denom¬ 
inators of the fractions which occur in Book 8. The resolutions used for B are (in h): 
1/360, i/eo, 1 / 30 , 1 / 24 , 1 / 15 , 1 / 12 , 1 / 10 , i/s, 1 / 7 , i/e, and i/fl For X i/eo and 1/7 were not 
use d, and 1 /20 an d 1/9 were added. The main resolution was set to Og = 1/12 h (follow¬ 
ing [Cunta, ll923|) and for a further run to Og = 1/24 h. The results for the portions of 
apparent resolutions are: 

- a. 


61%, a < ttg 
51%, a < Og 


- a. 


= 1/12 h: 

VL: a = a, 

'A: a = a, 

= 1/24 h: 

B: a = Og: 61%, a < a. 
X: a = Og: 55%, a < a., 


34%, a > Og: 4%, no fitting rounding: 1%; 
36%, a > Og: 9%, no fitting rounding: 4%; 

24%, a > Og: 14%, no fitting rounding: 1%; 
23%, a > Og: 19%, no fitting rounding: 4%. 


^The denominator 360 does not really occur but was introduce d because of the occurrence of the 
sum l/s -I- l/so = (see lStiickelberger and GraBhofil . I 2 OO 6 I . p. 883, n. 65). 
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The result indicates that Ug = V24h is the main resolution, but also Og = ^/i 2 \v comes 
into consideration. The result of Procedure ROUND applied with Og = 1/24 h shows that 
ca. 80% of the localities have only small absolute differences < 0.5 • 1/24 h = 1/48 h = 
1.25 min (cf. Eq. ([3|)) between Ai and A(Ai) converted by Eq. ([T]). Further, almost all 
Ai can be explained by conversion based on either Aa or Aa and a rounding. Hence, 
the vast majority of the differences between the given Ai and the converted A{Ki) are 
explicable by the precision of the Ai. Consequently, the location catalogue and Book 
8 show a high level of interdependence. 

The determined frequencies of Aa and Aa are: 

~ H: Aa: 107 (31%), Aa: 209 (61%), Aa or Aa: 25 (7%), no fitting rounding: 3 (1%); 
“ X: Aa 137 (62%), Aa: 63 (29%), Aa or Aa: 11 (5%), no fitting rounding: 9 (4%). 
This again reveals that Aa dominates in X and Aa in Vt (with a ratio of ca. 2:1). 

Figures [3] and 0] show the result of the classification for each locality in the Ptolemaic 
coordinate system (Aa: blue/circle, Aa: yellow/circle, Aa or Aa: green/square, no 
fitting rounding: red/triangle). In Q the dominant Aa occurs nearly exclusively in the 
east, and in the west also Aa occurs often. There are in fl 13 chapters in which Aa 
dominates: Europe 1, 3, 5, 6, 9, Asia 3, 5-11; in 2 chapters Aa dominates: Europe 7, 
Africa 2. In X Aa dominates in 8 chapters: Europe 2, 5, 7, 10, Africa 2-4, Asia 1; Aa 
dominates in 2 chapters: Europe 1, Asia 2. 

2.3 Comparison of H and X 

Finally, the Ai of 11 and X were compared. 343 values are equal and 143 values have 
differences, thereof 16 localities also have differences in A in the location catalogue 
so that 143 — 16 = 127 Ai remain with actual differences between 12 and X (37%). 
Besides the expected systematic difference of AAa = 30' = 2 min also other differences 
occur comparably often: 2.5 min = 1/24 h, 3 min = 1/20 h, 4 min = V^sh, and 5 min = 
1/12 h. This is explicable by the precision of the Ai. As an example Clunia (8.4.5) 
is considered: In H A is —3^15™, that is A = 11°00' converted with Aa = 60° to 
A = —3^16™ and rounded to the nearest 1/12 h. In X A is —3*'20™, that is A converted 
with Aa = 60°30' to A = —3^18™ and rounded to the nearest 1/12 h. The systematic 
difference of AAa = 2 min becomes 5 min here. 

Figure [S] shows the localities without (green/circle) and with differences in A be¬ 
tween H and X (Ptolemaic system). Those with differences are divided into localities 
with difference in A (yellow/triangle), without difference in A (red/triangle), and the 
remaining localities (magenta/square) which have no A in the location catalogue of X. 
Localities with differences disperse across the whole area; however, they often occur in 
local groups. In the west, there is a larger amount of differences. 

3 Length of the longest day 

The length of the longest day M increases from 12 h at the equator (4) = 0) to 24 h 
for 4) = 90° — e. Ptolemy derives the computation of 4) from M and vice versa in his 
Mathematike Syntaxis (MS) 2.3. The modern formulation of that relation between M 
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and d) is 


cos I -M 157h 


= — tan $ tan e 


( 4 ) 


(M in hours) with the obliquity of the ecliptic eH. In MS 2.6 Ptolemy gives a compi¬ 
lation of parallels with specific M including their latitudes (see Sect. Id.5|) . 

Ptolemy does not explicitly mention the e underlying the conversion between <i> and 
M in the Geography. In GH 7.6.7 he gives the ratio e : ~ that is the roughly 

rounded value 90° • Yis = 24° =: for e. In MS 1.12 Ptolemy states that the arc 

between the solstice points (i.e. 2e) i s approxirnately (11/83) • 360°. That is 47°42'40" 
so that £ becomes 23°51'20" =: £ni ( Manitiua . ll912L p. 44). Other appearances of £ni 
are MS 1.14 and 2.4, where Ptolemy uses this value for exemplary computations, as 
well as the “Table of Inclination” in MS 1.15, where occurs as the size of the arc of 
meridian with an ecliptic longitude of 90° (see Sect. 11. The correct value of Ptolemy’s 
ti me is 23°41' =: Sc so that £m has an error of about 11'. 

Rawlins! 1 1985h gives two groups of cities of Book 8, on the one hand 10 cities, whose 
M and d* could be related by Eq. through 23°50' ~ Em, on the other hand 7 cities. 


whose M and 4) could be re lated through 23°55' =: £h. Rawlins! attributes the latter 
value to Hipparchus (see also lRawlinsl . ll98dl so that the original conversion would have 
been done with that value before Ptolem y. A further hint to = 24° used in Ptolemy’s 
data was detected by Neueebau^ ( 19751 pp. 236, 245 ff.), who found £r underlying the 
data in the tables for the first and last visibility of the planets in MS 13.10. The 
question of the e underlying the Mi in Book 8 is considered in the following. 


3.1 Conversion of coordinates and accuracy 

A conversion between 4> and M was carried out for both directions M i—$ and 4) i—>■ M 
by means of Eq. (|T]), and the differences A4)i = 4>i — $(Mi) and AM^ = Mi — M(4>i) 
to the given values 4>i and Mi were computed. For M 4> with £ = £^ Fig. O shows 
the frequencies of differences A4?i lying within specified intervals (interval width 10'). 
The distribution has a bell-shaped form with its maximum at ~ 0. The majority of 
the absolute values of the differences is not larger than 1 / 2 °. 

As a measure of the consistency between given coordinates yi and converted coor¬ 
dinates y{xi) the median absolute deviation 


MAD = med(|?/i - y{xi)\) 
(i = 1.. .n) and the mean deviation 


( 5 ) 


1 

-^\y,-y{x^)\ 


( 6 ) 


are used here ( Hampel . 1974 : SachsL 19921 p. 337). Unlike the MD, the MAD is resistent 


^Trigonometric problems could be solved accurately in antiquity by using tabulations of chords as 
it c an be found in MS 1.11. On ancient trigonometry used by Ptolemy and described in his MS 
see iNeugebauerl lll97f]l . p. 21 ff.). 
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against outliers (i.e. abnormally large differences yi — y{xi) here). 

For both directions of conversion four conversions were carried out with 6^, Eh, £r, 
as well as experimentally with Sc- Table [T] shows the MAD, MD, and the number of 
coordinates n. Because of its resistent property the MAD is lower than the MD. In 
most cases value Em leads to the lowest measures of precision; only in the case of 
M $ the MAD based on Sh is I' (X) and 2' (SI) smaller than the MAD based on 
Em- Regarding $ M as the main direction, the results does not argue against Em 
underlying the data. 

A conversion of coordinates and a comparison with given coordinate values was 
also carried out for the Mi of Book 8 and the latitudes of the Canon. For e = Em 
the resulting MAD slightly increases in comparison to the location catalogue by 1' or 
0.1 min, respectively, and the MD increases noticeably by 6' or 0.4 min, respectively (see 
Table [T|) . Consequentl y the close relation between the M -data of X and the latitudes 
of the Canon stated by Stiickelberger and Grafiho^ ( 20061 p. 44) cannot be supported 
by means of the present computations. 

The size of the differences — $(Mi) can be explained by a rough resolution 

of the M-data and a disadvantageous propag ation of their errors (assuming that the 
original direction of conversion was 4) i—)■ M). ICunt d ( 19231 p. 104) assumes that the 
original resolution was 1/12 h. The largest denominators occurring in D among the 
fractions of the Mi are 150 {Palura (8.26.6)) and 60 (e.g. Ptolemais Herme (8.15.13), 
Omanon (8.22.12), Saue (8.22.15)); however, they are exceptional values. In X the 
largest denominator is 30. In both D and X the denominator 12 occurs often so that 
the main resolution of M-values seems to be 1/12 h = 5 min. From fractions like n/6 = 
2n/12 and n/4 = 3n/12 the precision is not known for sure, rougher resolutions than 
1/12 h are also possible. For 30° < 4> < 40° Fig. 0 exemplarily shows the data points 
($i,Mi) of D and X as well as the curve of the trigonometric conversion function 
M($) (derived from Eq. ([ 1 ]), e = Em)- The data points are arranged in horizontal 
groups, within which the M-values of different are equal. Most of these groups 
have a distance of 1/12 h so that their resolution is apparently 1/12 h (doted lines). Also 
Mi with other resolutions become apparent. The plot shows, too, that equal were 
converted to different Mi. That can be explained by a rougher precision, i.e. a = i/eh. 
Regarding only the lines of a = i/eh (solid), no with different Mi occur. 

The error of an M affects the corresponding latitude $(M) according to the value 
of that latitude. In general, the estimated effect of an error in x on the value of 
the function y{x) is 

Ay 




( 7 ) 


with the derivative dy/dx of y{x) with respect to x. Applied to Eq. (0) the variation 
A$ of 4> according to the error Am of M can be determined. 

Figure [5] shows the propagation of errors for the range of M in Book 8 (12 h < M < 
20 h) and for the resolutions a = i/soh, 1/30 h, and 1/12 h; the variations were set to 
Am = a/2 (Eq. ([3])). Most of the localities in Book 8 have an M < 17 h. For them 
a resolution of 1/12 h can cause a significant loss of accuracy, the variation A$ of 4) is 
> 10'. A<i, increases for M —>■ 12h to 2 / 3 °. Accordingly, an ancient conversion 4) 1 —>■ M 
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with a precision of 1/12 h (or more imprecise) caused a significant loss of accuracy 
in the latitude expressed by M. Even exceptionally high precisions as i/eoh could 
have a noticeable effect as the example of Ptolemais Herme shows: M = 1343/6oh, 
Am = 0.5 • i/eoh, A$ = 7'. Therefore, a large denominator in the fraction of M does 
not generally mean that the value is given with an unrealistically high precision. 


3.2 Distribution of the coordinate differences 


For the coordinate differences AM^ = Mi — M{^i) after the conversi on <i> ^ M w ith 
e = Em it was tested by means of the y^-test of goodness of fit (e.g. ISac y. I 1992 L p. 
420 ff.) whether they are normally distributed. The test statistic is in general 


k 

= '^{ni-npi)'^/{npi) (8) 

Z=1 

with the class number /c, the observed frequencies Ui of class i, the hypothetical prob¬ 
abilities Pi, and n = Null hypothesis of the test is the normal distribution 

assumption. Under this hypothesis follows a y^-distribution with f = k — 1 — u 
degrees of freedom, where u is the number of the unknown parameters. The two pa¬ 
rameters of the normal distribution, the location parameter p and the scale parameter 
a, were estimated by the arithmetical mean of the AM^ and by the root mean square 
deviation (with respect to the mean). A few outlying, very large differences, which are 
not characteristic for the distribution form, were excluded. According to the assumed 
crudest resolution a = 1/4 h of the Mi (see Sect. IXil) the limit for |AMi| was set to 
a/2 = 7.5 min (Eq. ([3])). The level of significance a used is 5%. The results are: 

- XI: n = 341, fc = 10, / = 7, /i = 0.06 min, a = 1.58 min, y^ = 26.9; 

- X: n = 216, fc = 10, / = 7, /i = —0.17min, a = 1.83 min, y^ = 27.5. 

The quantile is in both cases qi-aj = 14.1 so that the normal distribution assumption 
has to be rejected. 

The y^-test was repeated for the Laplace distribution. Its two parameters p and a 
were estimated by the median and by the MD (Eq. ([3]), but here with respect to the 
computed median). The results are: 

- XI: n = 341, k = 10, f = 7, p = 0.06 min, a = 1.16 min, y^ = 12.4; 

- X: n = 216, k = 10, f = 7, p = 0.03 min, a = 1.35 min, y^ = 7.6. 

For Xl and X y^ goes below the quantile qi-aj = 14.1. Accordingly it can be assumed 
that the differences AMi follow (at least approximately) a Laplace distribution. 


3.3 On the original direction of conversion 

Ptolemy states in GH 8.2.1 that the M-data are computed from latitudes. However, it 
i s known from a few l ocalit ies that Ptolemy adopted their A/-values from Hipparchus 
( Berggren and JonesLl200^ p. 29). Hipparchus attributed them to parallels with par¬ 
ticular lengths of the longest day [Meroe (8.16.9), Ptolemais Theron (8.16.10), Soene 
(8.15.15), Rhodos (8.17.21), and Massalia (8.5.7)). From these localities can be ex¬ 
pected that their latitudes in the location catalogue are derived from M in Book 8. 
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Rawlins (Il985l mentions the possibility of a conversion M i—$ for major cities; 

RawlinsI dlQSfJl Samples are given. Exemplarily, the localities shall be considered 


which iRawlina relates to Eh- The contrary conversion d* i—>■ M was done by means of 
Eq. dH) with Em; the result is shown in Tabled As can be seen, each computed M{^) 
agrees satisfactorily with the value in Book 8. Hence, the coordinates of these localities 
also support the direction of conversion d* i—)■ M. 

To provide evidence for the original direction of conversion it was examined how 
satisfactorily the coordinate differences and AMi = Mi — M{^i) 

after both conversions M $ and $ i—>■ M can be explained by the precision of the 
coordinates. Unfortunately, the precision of $ and M is not constant and not known 
in each case; however, a maximal resolution can be given for the <i>-values and a main 
resolution for the M-values. Thereby limits Ta$ and Tam for the differences and 
AMi can be given so that the numbers of differences was determined, which go below 
these limits. 

For the <i>-values of the poleis episemoi a maximal resolution of 1 / 2 ° can be as sumed , 
defin itely the actual resolution is often more precise (e.g. 1/12° or 1/6°, (see iMarxI . 
l201lll 'l. According to Eq. ([21) the limit Ta<i> is then 0.5 • 1/2° = 1/4°. For the M-values 
a main resolution of 1/12 h (or more imprecise; see Sect. IXTI) can be assumed so that 
Tam = 0.5 • 1/12 h = 1/24 h = 2.5 min was used as a limit. 

For both conversions M $ and $ 1 —>■ M the parameters Em and Eh were used. The 
result is: 


- Em: H: 69% |A$| < Ta-k, 88% |AM| < Tam; X: 66 % |A$| < La<s>, 82% |AM| < 
Lam] 

- Eh: H: 68 % |A$| < Ta<i„ 87% |AM| < Lam] X: 66 % |A$| < Ta-c, 81% |AM| < 
Lam- 

The test conditions for the conversions M $ and $ 1 —>■ M were not the same, for 
A<1> the limit was advantageously large and for AM small. Nevertheless, the amount 
of AMi going below their limit is noticeably larger than in the case of A<i>i. Thus, the 
results argue for a main direction of conversion of $ M. The resulting Mi of this 
conversion were round values, whose inaccuracies affect the A$i computed here after 
the back conversion M 1 —>■ d*. 


3.4 The obliquity of the ecliptic 

The question of which value e underlies the M-data was investigated more in detail 
chapter by chapter, assuming that the majority of the localities of a chapter were 
converted with the same value of e. An approach to determine the underlying e is 
an adjustment of e, which also leads to information about the accuracy of the result. 
Solving Eq. dH) for M gives the observation equation of the adjustment model 


Afi-b Ui =/(e) = 2^^ arccos (—tan$i tanE) , i = (9) 

11 ) 


with the observation Mi and its correction Vi. The solution of this non-linear ad- 
justment problem can be obtained iteratively by means of a linearisation (see e.g. 
Jag-er et al.l . 120051 p. 116). In Sect. 13.21 the y^-test of goodness of fit showed that the 
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differences AM^ with respect to a conversion with e^n follow rather a Laplace distribu¬ 
tion than a normal distribution. It shall be assumed here that the Laplace distribution 
can also be expected for errors of Mi whose conversion is based on e ^ Em. The max¬ 
imum likelihood estimation f or the Laplace dis tribution is the Li-norm estimation, 
which is applied here (see e.g. .Tager et H] . l2005l p. 117ff.). 

To estimate the accuracy of the adjusted parameter e, the term q = (A^A)“^ 
is needed, whose one-column design matrix A of the linearised adjustment problem 
contains the derivatives dfi{e)/de of fi with respect to e. Then the (1 — a) 100% 
confidence i nterval for e can be estimated according to iDielman and Pfaffenberger 
(119821 llQQflh by 

e ± Ae = e ± Za/ 2 ^\/q (10) 

with the quantile Za /2 of the standard normal distribution as well as with 


A = 


ivin 


( 11 ) 


Value m is the index of the median of the computed corrections Vi in the ordered 
sequence Vi. Values Vm+fj. and Vm-v are the corrections having a distance iz from Vm 
in the ordered sequence. 

The investigations in Sects. IXT] and l3.d] show that there are inaccurate Mi with larger 
errors, possibly partly corruplQ. The Li-norm estimation of e and the estimation 113 
of Ae have the advantage to be resistent against outlying observations. 

Due to their small number of observations (< 3) the chapters Europe 4 and Asia 12 
were excluded from the parameter estimation. A was computed using all Vi of the 
remaining chapters, v must to be small but larger than the half of the number u of 
adjustments of e; it was set to 11 -I- u/2 = 18 for X and to 18 -I- u/2 = 30 for fi. For 
the confidence interval of e the probability 95% was chosen. In summary, the results 
are 

- M. n = 338, M = 24, A = 0.7', smallest e: 23°46' ± 8' in Europe 10, largest e: 
24°07'±30' in Asia 11; 

- X: n = 217, u = 16, A = 1.0', smallest e: 23°30' ± 44' in Africa 4 (anomalously 
inaccurate) and 23°45' ± 11' in Europe 10, largest s: 23°56' ± 11' in Europe 9. 

With the exception of Asia 11 in D and Africa 4 in X the adjusted values of e deviate 
from Em in the order of some arc minutes (< 10'). Only in the case of Asia 3 in D 
Em lies outside of the estimated 95%-confidence interval (e = 24°02' ± 9'). Apart from 
this exception, the adjustment does not result in significant deviations from Em and 
accordingly the results do not argue against Em underlying the M-values in general. 

Moreover, Em need not be rejected as long as the differences AMi = Mi — M{^i) 
based on Em can be explained by rounding. Exemplarily Asia 3 in D with the significant 


'^E xamples for corrupt or roughly given M-values are Nisaia (8.23.6) and Arbis (8.25.14) in 
fStiickelberger and Mittenhubcif . 120061 . p. 205, 209) 

®In an Li-norm estimation of u unknown parameters, (at least) u corrections become zero. Th e choic e 
u' = 11 = I/ — tt/2 for X is reconcilable with the values used in iDielman and Pfaffenberged lll990h . 
According to the ratio (nfj — u^)/{nx — u-x) = l'6i was set to = l.6u' + u/2 = 18 -I- 12 = 30 
for Q. 
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deviation from e^n is considered. Each of its ten Mi can be explained by a kind of 
rounding as it can be found among the Mi of Book 8 (see Sect. I8.5|) . There are eight 
Mi explicable by a rounding to the nearest 1/12 h. Further, in the case of Thospia the 
precision could be 1/24 h (M($ = 39°50') = 14*'53.2™ si 14 21/24 h = 14 7/8 h = M) and 
in the case of Albana Yeh {M{^ = 45°50') = 15^36.7™ ~ ISYeh = 14 2/3 h = M). 
Apart from that, M of Artemita is not explicable by a conversion with Ej. (which is 
to be expected in view of the adjustment) in conjunction with any kind of rounding 
(M($ = 40°20') = 14^^57.7“, but M = 14*^55.0“). 

Finally, all Mi were used for an Li-norm-adjustment of one parameter e. For 11 is 
obtained e = 23° 52' ± 2' and for X e = 23° 52' ± 3'. Among the different values for e 
introduced here, only Em lies within both 95% confidence intervals. 


3.5 Conversion by means of linear interpolation 


It is much more probable that a conversion of coordinates was originally done by 
means of a tabulation an d not b y mean s of a single, trigonometric computation from 
Eq. dH) for each locality. Cunt j ( 19231 p. 101 ff.) assumes an underly ing tab ulation 
(Mj,$(Mj)) with a step size of Yi 2 h = 5min. Using this tabulation, Cunt d cannot 
explain each M-value in the examined chapters (8.3-8). An e xample is Gesoriacum 
(8.5.6) with 4) = 53°30' and M = IG^'SO™. The nearest 4? in Cuntz’s tabulation is 
53°39' which leads to M = 16''55™. Thus, Cuntd supposes that $ = 53°30' is not the 
original value underlying M. 

Another approach to explain the M-data in Book 8 is taking recourse to known 4>- 
and M-data given by Ptolemy. There are two datasets with latitudes and correspond¬ 
ing M-data in the Geography and the MS, these are the compilations of parallels in 
GH 1.23 and MS 2.6. It is conceivable that such data wer e used in conjunction with 
an interpolation, a method applied by Ptolemy (IPedersenl . 1201 iL p. 84). Considering 


the mentioned example of Gesoriacum^ Af(4>) computed from Eq. ([?]) is 16’°54™. The 
same result is obtained by means of the compilation of the parallels in GH 1.23 and a 
linear interpolation. The slope a of the 17th and 18th parallels in GH 1.23 is 


a — (Mj+i - Mj)/(4>j+i - 4)j) 

= (17‘'00“ - 16^^30“)/(54°00' - 51°30') = 0.2 V° • 


( 12 ) 


The approximated M of Gesoriacum is then 

M = 16^^30“ -b a • (^ - 51°30') = 16‘'54“ . (13) 

This value rounded to the nearest Ye is 16’'50™, the value in the location catalogue. 
As the plot in Fig. [7]shows, a precision of Yeh can be assumed for part of the M-data 
(see Sect. 13.11) . 

The accuracy of a linear interpolation of the data in GH 1.23 and MS 2.6 was deter¬ 
mined by computing the difference AM(4>) = M(4>) — M(4)) between the interpolated 
value M(4)) and the trigonometric computed value M(4>) (Eq. (H])). Figure [5] shows 
AM(4>) and additionally the differences AM^ = Mi — M(4)i) between the Mi of Book 
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8 and the M(^i) computed from the of the location catalogue (a few data points 
with I AM I > 7 min are not plotted). For d* > 54° the curves of the interpolated GH- 
and MS-data are equaj^. For d* < 54° the interpolation of the MS-data is very accurate 
(|AM| < 0.5 min); for $ > 54° the majority of the differences AM^ coincide with the 
curve AM($). This indicates that the Mi of Book 8 are based on an interpolation of 
the data in MS 2.6. The interpolation of the GH-data produces errors up to 3 min. 
Its curve AM($) does not coincide with the distribution of the differences AM^ in 
general. However, some data points lie next to the curve, and three data points coin¬ 
cide with the curve so that an interpolation of the data in GH 1.23 can be assumed 
at least in these cases (11: Harmuza (8.22.21), Aspithra (8.27.11); H and X: Ambrodax 
(8.21.17)). 

Finally, the trigonometric computation (e = Sm) and the linear interpolation were 
compared with regard to the amount of Mi explicable by their precision. For this 
purpose. Procedure ROUND was applied (see Sect. HI). As it was shown in Sect. o 
different resolutions come into consideration. The possible resolutions of M were chosen 
according to the denominators of the fractions occurring in Book 8. In addition to the 
standard resolution Og = Vi 2 h there are in 11 the more precise a (in h) i/is, 1 / 20 , 1 / 24 , 
1 / 30 , i/eo, and 1/150 as well as the less precise 1 / 10 , i/s, i/e, and 1 / 4 . The more precise 
a in X are (in h) 1/15 and 1/30 as well as the less precise i/s, i/e, and 1 / 4 . The result is 
given in Table [31 The more precise the resulting resolutions are, the more probable the 
conversion can be considered. The trigonometric conversion and the interpolation of 
the MS-data have nearly the same result; the interpolation of the data in GH 1.23 has 
a worse result than the other two methods. Thus, a conversion by means of the data 
in MS 2.6 is very probable. The amount of Mi being explainable with a resolution of 
1/12 h is 72% for H and 69% for X. 

The results in Table [3] show that in the case of 11 84% and in the case of X 72% of 
the localities have only small absolute differences < 0.5 • 1/12 h = 1/24 h = 2.5 min (cf. 
Eq. ([31)) between Mi and M{^i) converted by means of the data of MS 2.6. Almost 
all Mi (ca. 95%) can be explained by a conversion (data of MS 2.6 or Eq. ([1])) and 
rounding. Hence, the majority of the differences between Book 8 and the location 
catalogue are explicable by the precision of the M-data. Accordingly, both datasets 
show a high level of interdependence. 

There is a lower percentage of Mi (3% in H, 6% in X) whose differences AM^ 
are inexplicable by one of the types of rounding used. These as well as lower precise 
explicable M-values occur in nearly all chapters of Book 8. Possible reasons for Mi 
and M{^i) not fitting are 

(1) M is rounded incorrectlj0 or is slightly erroneous due to a propagation of rounding 
errors in intermediate results; for example possibly Asturica August (8.4.5) with 
M = 15‘'25“ (in H and X) and M($) = 15‘'22.2“ si 15*'20“; 


®There are more parallels in MS 2.6 than in GH 1.23 from 45° latitude, but from 54° to 63° there 
are some interpolated values in MS 2.6 so that a linear interpolation of the data in MS 2.6 and in 
GH 1.23 leads to the same result for latitudes > 54°. 

known inconsistence of this kind, but in latitude, occurs in GH 1.23 in the latitude 4> = 43°05^ 
of the 13th parallel. Since M of this parallel is 15^15™, the correct value for $ should be 43° 00^ ~ 
4>(M) = 43°01^ which was given by Ptolemy in MS 2.6. 
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(2) M is wrong due to a corruption; for example possibly Palmyra (8.20.10) in with 
M = 14^40“ and M($) = 14*'19“ Ri 14^20“; 

(3) 4) is wrong due to a corruption; for example Osika (8.19.9) in X with $ = 45°30', 
M = 15‘’27.5“, M($) = 15‘'34“ (in contrast fl: 4> = 44°45' and Canon: 4> = 
44° 30'). 

It was tested for the remaining inexplicable Mi, whether a conversion with Sh or 
Er brings an explanation of their fractions. For this purpose Procedure ROUND was 
repeated. As a result, a few Mi can be explained by a rounding, but there still remain 
inexplicable values. However, it is not necessary to assume another value than e^; 
each of the problematic Mi can be explained by Error 1: Asturica August: see above; 
Luppia (8.6.3): M($) = 16*'44™ became 16^50“ instead of 16^40™ (resolution V^h); 
Ravenna (8.8.5): as Asturica August; Capua (8.8.6): M(4)) = 15''02™ became 15*'05™ 
instead of 15*'00“; Hallos (8.17.44): M(4>) = 14*'32.8“ became 14^34“ instead of 
14^32™ (resolution i/soh); Medaba (8.20.20): M($) = 14^02™ became 14*'05™ instead 
of 14‘'00“. 

4 Distance from the summer solstice 

At places located between the tropics the sun is at the zenith two (for |4)| < e) or one 
(for 14)I = e) times within the year. For this date, the distance S of the sun from the 
summer solstice on the ecliptic is given in Book 8. The quantity I = 90° — S' is the 
ecliptic longitude, reckoned from the vernal equinox. There are about 60 localities in 
Book 8 with given S. In MS 1.14, Ptolemy explains the computation of the distance of 
the sun from the equator (declination) by means of H. When the sun is at the zenith, 
its declination equals the latitude 4) of the locality. Thus, the calculational relations 
can be used to determine I or S from 4>. The relation between S and the latitude 4) is 

sin4> = sine cosS . (14) 

In MS 1.15, Ptolemy gives a list of tabulated ^{Ij) for Ij = 1°,...,90° in steps 
of 1° (“Table of Inclination”). A further source of S-data is MS 2.6, where Ptolemy 
arranges data about specific parallels (see Sect. I3.5|l : however, it includes only 5 data 
pairs {^j,Sj). For the parallels with 4) < e S' is given there. In the following both 
compilations of S-data in the MS are taken into consideration besides Eq. as a 
possible basis of the S-data in Book 8 of the Geography. 

4.1 Conversion of coordinates and accuracy 

Firstly, the precision of the Si is considered to get an idea, which size of coordinate 
differences are to be expected after a conversion between S and 4). The largest de¬ 
nominator of the fractions of the Si is 6 in 11 and 8 in X; the main resolutions seem 
to be 1 / 3 ° and 1/4°- Solving Eq. (fill) for 4>, the variation A$ of 4) due to an error As 
in S can be estimated by Eq. ©• The computation of propagation of errors was done 


®In MS 1.14, the ecliptic longitude I is reckoned from the autumnal equinox. 
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exemplarily for the resolutions 1/6° and 1/2°. The variations were set to the expected 
maximal rounding errors (Eq. @): A 5 = 0.5 • 1/6° = 1/12° and As = 0.5 • 1/2° = 1/4°. 
Figure [S] shows the result. The effect of a variation in S' on $ is very low (in comparison 
to M and $). For an inaccurate resolution of 1/2° the expected maximal error in $ 
increases only to 6'; it goes to zero for S —>■ 0 (i.e. <&—)>£). 

The conversion between S and $ was carried out for both directions S 1 —$ and 
$ S by means of Eq. m with the four different e introduced in Sect. [31 and 
the differences — $(Si) and ASi = Si — S($i) to the given values <1)^ and 

Si were computed. Since X ends in GH 5.13, only 16 S-values can be considered for 
the comparison of Book 8 and the location catalogue (these are Syene (8.15.15) and 
Chapter 16 (Asia 4)). From the differences between the converted and the original 
coordinates, the MAD and the MD were computed by Eqs. (O and (jS]); the results are 
given in Table [ffl. 

In the case of X, the MAD = 2' for $ S(£i„) is very small and goes below the 

expected result with regard to the assumed main precisions of 1/3° and 1/4°. In contrast, 
the MD = 31' is considerably larger. It is increased by four large differenced^. The 
MAD for £ £m is significantly larger than for £„; that argues for underlying the 
5' -data of X. _ 

Stiickelberger and GraBhi^ ( 20061 . p. 44) mention that the ^-data of X is consid¬ 


erably more precise that the .S-data of D. Using a conversion based on Eq. (11411 this 
can be confirmed. In the case of U, the MAD for $ 1 —>■ 5'(£in) is 72'. There are many 
differences ASi > 1° which cannot be explained by the resolution of the Si. A wrong 
value for e underlying the reconversion can be excluded; the other (eligible) values 
do not lead to considerably smaller coordinate differences. Thus, either corrupt data 
must be assumed for U or that the conversion done here is wrong. Regarding the known 
sources of 5'-data, the tabulation in MS 1.15 has a step size of 1° for I and Pto l emy’s 


Handy Tables include a similar tabulation with a step size of 3° ( Neugebauerl . Il97,^ 


p. 979). Thus, the usage of these tabulations (without interpolation) would result in 
S'-values with fractions equal to zero. That does not hold for the data in Book 8. In 
Fig. (inithe differences ASi = Si — 5'($i) are plotted against the latitudj^. The plot 
indicates a dependence between AS and 4) of H; regionally AS increases or decreases 
with increasing 4). This suggests that the 5'-data in H are not determined by means of 
the conversion dD. 

Table H] also gives the MAD and the MD for the differences A4)i and ASi computed 
from the 4>i and converted 5'(4'i) of the Ganon. In the case of fl the results do not 
change significantly in comparison to the location catalogue of the Geography. In the 
case of X the Ganon enables a comparison for all ^-data (in contrast to the incomplete 
location catalogue of X). Using the Canon, the small MAD = 3.5' for 4> 5'(£ni) 

confirms the high accuracy of the S'-data which has already been observed for the 


®In the case of <J> 1 —t S{£c) the number n of the converted latitudes is smaller because values $ > £c 
cannot be converted by Eq. imi. Arbis (8.25.14) is excluded from the computations because of its 
exceptionally wrong S = 0° in Q. 

^^These localities are Syene (8.15.15) and Autolalai (8.16.3) with the crude value 5 = 0° for $ = 
23°50^ as well as Dere (8.16.12) and Mosylon (8.16.13). 

Arbis (8.25.14) with very large |A5| is excluded in the case of fl. 
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small sample size of the location catalogue. The large MD of 50' shows, as in case of 
the location catalogue, that there are also some latitudes which do not agree with the 
respective S. 

Due to the incomplete location catalogue of X the missing latitudes were experi¬ 
mentally substituted by the values of D as well as of the Canon to compute the 
In the case of D there are eight localities more than in the case of the Canon with 
large differences |A5'i| > 1°. Figure [Till shows the ASi based on the latitudes of 
the Canon. In the case of X there is no comparable dependence between AS and $ 
as detected for D. Most of the differences are near zero. Consequently, the coordinate 
differences of X reveal a close connection between its S'-data and the latitudes of the 
Canon. 

Furthermore, Fig. [ITll shows the curves of the differences A,S(<I>) = 5(<i)) — ^(d)) be¬ 
tween the trigonometrically computed 5'(d>) (Eq. (flUl l and the value ^(d)) determined 
by linear interpolations: 1) of the tabulation in MS 1.15, 2) of the tabulation in the 
Handy Tables and 3) of the data of the parallels in MS 2.6. The interpolations of 
the data in MS 1.15 and the Handy Tables lead to accurate results (|A,S| < 2.3' and 
|A,S| < 0.5° respectively). Their error curves run close to zero so that they coincide 
with the majority of the differences ASi occurring in X. Consequently, a conversion 
by means of one of both tabulations can be assumed for the predominant amount of 
accurate Si in X. Because of the small number of data pairs (d)j, Sj) in MS 2.6, the 
interpolation of these data is significantly erroneous, especially for d> > 20°. However, 
the curve AS fits the differences ASi of some localities in H for d> > 16° well. This 
reveals the possibility of another type of interpolation underlying the data of H, see 
Sect. 14.21 

For the S'-data of X and the auxiliary latitudes of the Canon, Procedure ROUND 
(see Sect. [a was applied to explain as many Si as possible by rounding. Because of 
the small MAD, the standard resolution was set to Os = = 7.5'. The less precise 

resolutions are a = 1/6°, 1/4°, 1 / 3 °, 1 / 2 °, and 1°. Asa result 42 of 59 Si can be explained 
by a rounding: a = as'. 15, a = 1 / 6 °: 20, a = 1/4°: 3, a = 1 / 2 °: 3, and a = 1°: 1. Among 
the remaining 17 Si there are 5 with the obviously approximative value S = OTu. 
Furthermore, 2 Si can be explained by a linear interpolation of the data in MS 2.6, the 
concerning localities are Argyre (8.27.10) and Kattigara (8.27.13). Their $ is 8°30', 
the interpolated S is 68°49', that is « 68 3/4°, the S in Book 8. The remaining 10 
problematic Si probably have corruptions; for example, Dere (8.16.12) has d* = 11°00' 
and S = 62 3/4°, but the converted value is 5'(<i>) = 61°51' ~ 613/4° so that an error 
of 1° in S is possible. 


l^Localities with |AS| > 1° occurring in ft but not in the Canon are: Kane (8.22.9), Germ (8.22.10), 
Sapphara (8.22.16), Chaberis (8.26.5), Palura (8.26.6), Patala (8.26.10), Barbarei (8.26.11), Tosalei 
(8.27.5), Triglyphon (8.27.7), and Aspithra (8.27.11). 

^^The tabulation of the Handy Tables is represented here by the data of the table in MS 1.15. To 
obtain a step size of 3° for I each third row was used. 

^^These localities are Syene (8.15.15), Autolalai (8.16.3), Gerra (8.22.10), Harmuza (8.22.21), and 
Tosalei (8.27.5). 
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4.2 Linearly converted data in fl 

A few localities in Book 8 of have accurate S'-values with IAS"! < 15', they are Geira 
(8.16.6), Garama (8.16.7), Napata (8.16.8), Muza (8.22.6), Karminna (8.22.23), and 
Barbarei (8.26.11). The localities Syene (8.15.15), Berenike (8.15.19), and Autolalai 
(8.16.3) have d* = 23°50' and the approximative S = 0°. Arbis (8.25.14) has the 
erroneous 5 = 0°. The 5-values of most of the other localities can be explained by a 
conversion of $ by means of a piecewise linear function. The data reveal three regions 
for which a linear conversion was done separately. These regions cannot be made out 

well in the plot of the differences ASi in Fig. (TUlbut by means of a plot of 5 against 

4), see Fig. [TT] Obviously, the limits between the three regions are ^12 ~ 12°30' and 
4)23 ~ 20°15'; they correspond to a division of the quadrant of the ecliptic into three 
equal parts with the range of « 30°. Hence, the piecewise linear function of conversion 
is 

f 5(($) = Old) - 1 - bi for 0° < d) < d>i 2 
5'(d)) = < 52 (d>) = a 2 d) - 1 - 62 for d)i 2 < d> < d )23 (15) 

[ 53 ( 4 *) = a 3 d)- 1-63 for d )23 < d> < e 

To define a linear functions two points are necessary. Two points (d>, 5) are known: 
(0°,90°) and (e, 0°). Moreover, very probably Badeo (8.22.4) with d>B = 20°15' and 
5 b = 30° was used as end point of 5((d)) and starting point of 52 (d>) as well as 
Pityndra (8.26.19) with $p = 12°30' and 5p = 60° was used as end point of 52 (d)) and 
starting point of 53 (d)). 5 b and 5p are inaccurate, but conspicuously round values, 
possibly intensionally approximative to simplify the computation^^. It is obtained 


with (0°,90°) and Pityndra: 

ai = (5p - 90°)/(d>p - 0°) Ri-2.4 (16) 

= 90° - oi • 0° = 90° . (17) 

It is obtained with Badeo and Pityndra: 

0-2 = (5b — 5p)/(d)B — d)p) R! —3.9 (18) 

62 = 5b - 02 • d-B = 108°58' R! 109° . (19) 

It is obtained with Badeo and (e, 0°): 

03 = (0° - 5B)/(e - $b) ~ -8.37 ~ -8.3 (20) 

63 = 0° - 03 • e = 197°49' r: 198° , (21) 


where e = 23°50' is used, the latitude of localities with 5 = 0° in Book 8 (see above). 

5'(d)) is plotted in Fig. [TT] together with the given data points (4)^,5^) and the 
trigonometric conversion function 5(d)) based on Eq. ([TTI) . S' fits the inaccurate data 
points very well. Especially S'^ is a bad approximation for the trigonometric function 


^®An accurate S for ‘I>b exists in O. Napata (8.16.8) has <I> = 20°15' and S = 31°10b the value 
computed from Eq. l|Ti|l is S(<h) = 31°09'. 
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S'(<i>) based on Eq. (fTH) . In Tabled] S'i, the computed as well as the remaining 

difference AS"' = Si — S'{^i) are given for the 47 localities which have not been 
discussed so far. The differences AS- are explicable in almost every case by a reasonable 
rounding to the nearest 1°, 1 / 2 °, 1 / 3 °, or 1 / 4 °. There remain 10 localities whose S"(<i)) 
cannot be rounded to the given S in this way. The reason could be peculiar roundings 
or scribal errors, also in an originally trigonometric converted S. For example, a linear 
conversion and a confusion of 1/3 = 7 and 2/3 = yo can be assumed for Barygaza 
{S'{^) = 41°24' Ki 41 1 / 3 ° ^ 412 / 3 ° = S) as well as for Maagrammon. A trigonometric 
conversion and a scribal error can be assumed for lazeitha, which has S = 47°40' in 11 
and the accurate value 48°40' in X. 


Appendix: Changes to and errors in the data source 


The following changes of coordinates values underlie the present investigations. Changes 
of known inconsistent values are: $ of Autolalai (4.6.24) in X was changed from 27°20' 
to 23°50' (as 12 and Canon). Argyre (8.27.10) and Kattigara (8.27.13) were positioned 
on the northern hemisphere instead of the southern one. (The positive latitudes equal 
the mo re plausible val ues in the Canon. On the negative latitudes of Chinese places 
see also Rawlinsl (120081 on . 47).) Further chan ges are: Palmyr a (8.20.10): M = 14^40™ 
(also iNobbel . Il843-184^ to M = 1 4*^20™ (Rawliiii . 2009h and Apameia (8.20.9): 
A = i/sh to 5/8h (jNobbel . Il843-1845 : Rawlins . 2009ll . 

As a support for further investigations of the coordi nates of Book 8, the errors 
shall b e given here, which were found in the data base of lStiickelberger and GraBhofl 
( 2 OO 6 II . Data which are correct in the printout but erroneous in the data base are given 
in Table O 
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Table 1: Median absolute deviation MAD (Eq. Q) and mean deviation MD (Eq. ([5])) 
after conversions of coordinates of the fl-recension, the X-manuscript, and 
the Canon; M i—>■ $: deviations — d)(Mi) between the latitudes 
computed from the lengths of the longest day Mi of Book 8, and the <1)^ of 
Books 2-7/the Canon; $ i—)■ M: deviations Mi — M{^i) between the M{^i), 
computed from the of Books 2-7/the Canon, and the Mi of Book 8; S' $ 
and d) I—)■ S: as M I—$ and $ M but for the distances Si of the sun from 
the summer solstice of Book 8; n is the number of coordinates; the values of the 
obliquity of the ecliptic underlying the conversions are: Em = 23°51'20", £h = 
23°55',£r = 24°, £c = 23°4T 


Conversion 

0 B.2-7/B.8 

Canon/D B.8 

X B.2-5/B.8 

Canon/X B.8 


e 

MAD 

MD 

MAD 

MD 

MAD 

MD 

MAD 

MD 

n 


344 

344 

346 

346 

220 

220 

354 

354 

M $ 

'^m 

9.7' 

13.0' 

10.7' 

17.5' 

9.6' 

13.6' 

10.7' 

19.8' 


Eh 

7.9' 

13.6' 

11 .6' 

18.1' 

8.5' 

14.4' 

10.7' 

20.3' 


£r 

11.4' 

15.9' 

12 .8' 

20 .1' 

12 .0' 

17.3' 

11 .8' 

22 .2' 


£c 

14.3' 

18.4' 

15.2' 

22.3' 

15.4' 

18.6' 

15.8' 

25.1' 

4) M 


0.9” 

1.3” 

1 .0” 

1.7” 

1 .0” 

1.5” 

1 .1” 

1.9” 


£h 

0.9” 

1.4” 

1 .2” 

1 .8” 

1 .1” 

1 .6” 

1 .1” 

2 .0” 



1 .2” 

1 .6” 

1.5” 

2 .1” 

1.5” 

2 .0” 

1.4” 

2 .2” 


£c 

1.5” 

1.9” 

1 .6” 

2.3” 

1 .8” 

2 .1” 

1.7” 

2.5” 

n 


58 

58 

59 

59 

13 

13 

58 

58 

S $ 


18.7' 

22 .0' 

19.3' 

27.3' 

0.5' 

3.5' 

l.T 

8 .6' 


£h 

21 .2' 

23.1' 

21.7' 

28.0' 

3.2' 

5.4' 

3.1' 

10 .2' 


£r 

22.4' 

24.8' 

22.4' 

29.4' 

6.9' 

8 .T 

6 .6' 

12 .6' 


£c 

13.1' 

21 .0' 

15.0' 

27.3' 

7.5' 

9.5' 

7.5' 

12.9' 

4) S 


72.4' 

104.6' 

86 .2' 

117.5' 

1.5' 

31.2' 

3.5' 

50.0' 


£h 

81.5' 

117.4' 

Csj 

00 

124.1' 

11 .6' 

57.T 

11.3' 

64.3' 


£r 

00 

131.7' 

79.6' 

138.6' 

23.5' 

82.2' 

21 .8' 

80.4' 

n 


55 

55 

54 

54 

11 

11 

54 

54 


£c 

57.4' 

90.T 

72.7' 

110.3' 

23.0' 

32.2' 

22 .2' 

57.T 
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Table 2: Examples of the conversion of the latitude d* of the location catalogue into the 
length of the longest day M($) by Eq. ([I]) with e = e^n', M is the coordinate 
value of Book 8 


Locality 

$ 

[1 

[h] 

M 

[h,min] 

M(d>) 

[h,min] 

Arbela (8.21.3) 

371/4 

145/8 

14,17.5 

14,17.2 

Athenai (8.12.18) 

371/4 

145/8 

14,17.5 

14,17.2 

Karchedon (8.14.5) 

322/3 

141/5 

14,12.0 

14,11.8 

Lysimachia (8.11.7) 

411/2 

15 1/12 

15,05.0 

15,04.3 

Nikaia (8.17.7) 

4111/12 

151/8 

15,07.5 

15,07.1 

Rhodos (8.17.21) 

36 

141/2 

14,30.0 

14,29.9 

Tarentum (8.8.4) 

40 

1411/12 

14,55.0 

14,54.2 


Table 3: Estimated percentages of different resolutions a of the lengths of the longest 
day M of Book 8 (fl-recension and X-manuscript) in the case that they orig¬ 
inate from a trigonometric conversion (Eq. ([5])) or from a linear interpolation 
of the M-data in MS 2.6 and GH 1.23; last row: M does not fit with $ 







X 


trig. 

interp. 

MS 

interp. 

GH 

trig. 

interp. 

MS 

interp. 

GH 

a < 1 / 12 h 

12 

12 

11 

3 

3 

2 

a = h 

72 

72 

68 

71 

69 

65 

a > 1 / 12 h 

12 

12 

18 

21 

22 

24 

M ^ M($) 

4 

3 

4 

5 

6 

9 
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Table 4: 0-recension: Linear conversion of the latitude d* of the location catalogue 
into the distance S"($) from the summer solstice by means of the piecewise 
functions and ^^(d)) (Eq. (fT^ l and remaining difference AS" = 

S — S'(d>) with respect to S of Book 8 


Locality 

d) 

[M 

S 

[M 

S' 

r/] 

AS' 

['] 

S((d>) 





Adulis (16.11) 

11,20 

62,00 

62,48 

-48 

Dere (16.12) 

11,00 

63,45 

63,36 

9 

Mosylon (16.13) 

9,00 

68,45 

68,24 

21 

Aromata (16.14) 

6,00 

76,00 

75,36 

24 

Okelis (22.7) 

12,00 

61,15 

61,12 

3 

Arabia (22.8) 

11,30 

62,20 

62,24 

-4 

Palura (26.6) 

11,30 

62,00 

62,24 

-24 

Takola (27.3) 

4,15 

80,00 

79,48 

12 

Zabai (27.4) 

4,45 

78,40 

78,36 

4 

Argyre (27.10) 

8,30 

70,00 

69,36 

24 

Kattigara (27.13) 

8,30 

70,00 

69,36 

24 

Nagadiba (28.3) 

8,30 

70,00 

69,36 

24 

Talakori (28.4) 

11,40 

62,00 

62,00 

0 

Maagrammon (28.5) 

7,20 

72,40 

72,24 

16 






larzeitha (16.4) 

15,30 

47,40 

48,33 

-53 

Thamondokana (16.5) 

17,00 

43,00 

42,42 

18 

Meroe (16.9) 

16,25 

45,00 

44,58 

2 

Ptolemais Theron (16.10) 

16,25 

45,00 

44,58 

2 

Town of Pudnos (22.5) 

16,30 

44,40 

44,39 

1 

Kane (22.9) 

12,30 

60,45 

60,15 

30 

Mara (22.11) 

18,20 

37,30 

37,30 

0 

Omanon (22.12) 

19,20 

32,00 

33,36 

-96 

Menambis (22.13) 

16,30 

45,00 

44,39 

21 

Sabbatha (22.14) 

16,30 

45,00 

44,39 

21 

Sabe (22.15) 

13,00 

58,00 

58,18 

-18 

Sapphara (22.16) 

14,30 

52,30 

52,27 

3 

Isle of Sarapias (22.18) 

17,30 

41,00 

40,45 

15 

Simylla (26.3) 

14,45 

51,20 

51,29 

-8 

Muziris (26.4) 

14,00 

54,30 

54,24 

6 

Chaberis (26.5) 

15,45 

47,36 

47,35 

1 

Barygaza (26.12) 

17,20 

41,40 

41,24 

16 

Ozene (26.13) 

20,00 

31,00 

31,00 

0 

Baithana (26.14) 

18,10 

38,15 

38,09 

6 

Hippokura (26.15) 

19,10 

34,20 

34,15 

5 
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Locality 

$ 

[“/] 

S 

[M 

S' 

r/] 

AS" 

['] 

Karura (26.16) 

16,20 

45,20 

45,18 

2 

Modura (26.17) 

16,20 

45,20 

45,18 

2 

Orthura (26.18) 

16,20 

45,20 

45,18 

2 

Triglyphon (27.7) 

18,00 

39,00 

38,48 

12 

Mareura (27.8) 

12,30 

60,00 

60,15 

-15 

Aspithra (27.11) 

16,15 

45,40 

45,38 

2 

Sinai (27.12) 

13,00 

58,00 

58,18 

-18 


Gerra (22.10) 

23,20 

4,20 

4,20 

0 

Harmuza (22.21) 

23,30 

3,00 

2,57 

3 

Samydake (22.22) 

22,40 

10,00 

9,52 

8 

Patala (26.10) 

21,00 

23,50 

23,42 

8 

Tosalei (27.5) 

23,20 

4,20 

4,20 

0 

Tugma (27.6) 

22,15 

13,00 

13,19 

-19 


Table 5: Errors in 


the data base of Stiickelberger and Grafihofi ( 2006ll 


ID 

Name 

n/x 

Error 

Correction 

8.3.7.1 

Eboracum 

X 

link to 2.3.16 

2.3.17 

8.5.7.3 

Narbo 

X 

M = 16^15“ 

15^15“ 

8.5.7.4 

Vienna 

D 

A = 21^30“ 

21^15“ 

8.5.7.4 

Vienna 

X 

A = 21^15“ 

21^30“ 

8.5.7.4 

Vienna 

X 

M = 17^15“ 

15^30“ 

8 .8.6.1 

Aquileia 

D, X 

M = 15^05“ 

15^30“ 

8.13.8.1 

lol Caesarea 

X 

M = 14*'20“ 

14*'15“ 

8.17.32.1 

Perge/Aksu 

X 

A = 0*'30“ 

0*'05“ 

8.17.46.1 

Adana 

D 

link to 5.8.8 

5.8.7 

8.23.12.1 

Maruka 

Q 

CO 

II 

3 5/6h 

8.25.7.1 

Ortospana 

X 

M = 15*'20“ 

14h2om 

8.27.7.1 

Trilingon 

n 

A = 7*'16“ 

6*'16“ 
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Figure 1: Overview of the chapters in Book 8 (E.: Europe, Af.: Africa, As.: Asia) 
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Figure 2: Frequency of differences AAi = A^ — A{Ai) between the longitudes A^ of the 
location catalogue and the longitudes A(Ai) computed from the time differ¬ 
ences Ai to Alexandria of Book 8 by means of Eq. ([T]) for the 0-recension 
(left) and the X-manuscript {right); Aa is the longitude of Alexandria un¬ 
derlying the conversion 
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Figure 3: Classification of the localities in Book 8 of the fl-recension according to the 
longitude Aa of Alexandria which underlies a conversion of the longitudes 
of the location catalogue into the time differences Ai to Alexandria of Book 
8; blue/circle: Aa = 60°30', yellow/circle: Aa = 60°00', green/square: Aa or 
Aa, red/triangle: A does not fit with A 


0 ” 30 ” 60 ” 90 " 



Figure 4: Classification of the localities in Book 8 of the X-manuscript according to the 
longitude Aa of Alexandria which underlies a conversion of the longitudes A^ 
of the location catalogue into the time differences Ai to Alexandria of Book 
8; blue/circle: Aa = 60°30', yellow/circle: Aa = 60°00', green/square: Aa or 
Aa, red/triangle: A does not fit with A 
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Figure 5: Differences in the time differences Ai to Alexandria between Book 8 of the 
f2-recension and the X-manuscript; green/cirlce: no difference in A] other: 
difference in A] red/ triangle: no difference in the longitude A of the location 
catalogue, yellow/triangle: difference in A, magenta/square: no comparison 
of A 
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X 




Figure 6: Frequency of differences = <1)^ — ^{Mi) between the latitudes of the 
location catalogue and the latitudes $(Mi) computed from the lengths of 
the longest day Mi of Book 8 by Eq. (|3|) for the D-recension {left) and the 
X-manuscript {right) 
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Figure 7: Plot of the lengths of the longest day Mi of Book 8 against the latitudes 

of the location catalogue for the fl-recension and the X-manuscript, trigono¬ 
metric conversion function M($) (derived from Eq. (|3|)); the horizontal lines 
show the position of M-values with resolutions of a = i/eh and a = 1/12 h 
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Figure 8: Propagation of errors for the computation of the latitude <1> from the length 
of the longest day M (left) and the distance S from the summer solstice 
(right) for different errors Am and As 
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Figure 9: Differences AM = M — M($) between the length of the longest day M(<1>), 
trigonometrically computed from the latitude (Eq. ([I])), and the M-data 
of Book 8 (ri-recension and X-manuscript) as well as M resulting from an 
linear interpolation of the M-data of GH 1.23 and MS 2.6 
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Figure 10: Differences AS” = S — 5'($) between the distances S'(d)) from the summer 
solstice, trigonometrically computed from the latitude $ (Eq. (fUll 'l. and 
the S'-data of Book 8 (D-recension and X-manuscript) as well as S resulting 
from an linear interpolation of the iF-data of MS 1.15, MS 2.6, and Handy 
Tables (AS' « 0 for MS 1.15 and Handy Tables) 
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Figure 11: H-recension: Plot of the distances Si from the summer solstice of Book 8 
against the latitudes <1)^ of the location catalogue, trigonometric conversion 
function 5'(<i)) (derived from Eq. (fHl) '). linear conversion functions S[, 
and S'^ (Eq. (fT^ l 
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